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Abstract 



We investigate Bloch electrons in two dimensions subject to constant 
electric and magnetic fields. The model that results from our pursuit is gov- 
erned by a finite difference equation with a quasienergy spectrum that inter- 
polates between a butterfly-like structure and a Stark ladder structure. These 
findings ensued from the use of electric and magnetic translation operators. 
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We consider the problem of an electron moving in a two-dimensional lattice in the pres- 
ence of applied electric and magnetic fields. We refer to this as the two-dimensional electric- 
magnetic Bloch problem (EMB). The corresponding magnetic Bloch system (MB) has a long 
and rich history. An important early contribution to the analysis of the symmetries of the 
MB problem was made by Zak |I] , who worked out the representation theory of the group of 
magnetic translations. The renowned Harper equation was derived assuming a tight-binding 
approximation M, and Rauh derived a dual Harper equation in the strong magnetic field 
limit. The studies of Hofstadter and others Q of the Harper equation spectrum have since 
created an unceasing interest in the problem because of the beautiful self-similar structure of 
the butterfly spectrum |J. A remarkable experimental realization of the Hofstadter butter- 
fly was recently achieved, not for an electron system, but in the transmission of microwaves 
through an array of scatters inside a wave guide ||. The symmetries of the EMB prob- 
lem were analyzed some time ago by Ashby and Miller 0, who constructed the group of 
electric-magnetic translation operators, and worked out their irreducible representations. In 
this paper we utilize the properties of the electric-magnetic operators in order to derive a 
finite difference equation that governs the dynamics of the EMB problem. The numeri- 
cal solution of this equation displays an interesting pattern which interpolates between a 
butterfly-like structure and a Stark ladder structure. 

For the purposes of our research we consider the motion of an electron in a two- 
dimensional periodic potential, subject to a uniform magnetic field B perpendicular to the 
plane and to a constant electric field E, lying on the plane according to E = E(cos6,sin8) 
with 9 the angle between E and the lattice x\— axis. The dynamics of the electron is governed 
by a time-dependent Schrodinger equation that for convenience is written as 

^(t,f) = 0, (1) 

where 7r = po + A and tt = p + A , with the momentum operator p^ = (ij^, —iV^j. Units 
have been chosen here in which h = c = e = 1. Where necessary we use a covariant 
notation with space-time three vectors x^ = (t,r); \i = 0, 1,2. Eq. ([!]) can be consider as 
an eigenvalue equation for the operator S with eigenvalue 0. We adopt a gauge-independent 
procedure, thus the gauge potential is written as: 

A = (P-l)r-E, 

A 1 = (a - 1/2) Bx 2 - pEit, 

A 2 = (a + l/2)B Xl - (3E 2 t. (2) 

This potential yields the correct background fields independent to parameters a and (3. A 
general potential can be represented by its Fourier decomposition, however for simplicity we 
shall consider the potential 

U (xi, x 2 ) = U\ cos (2nxi/a) + U 2 cos {2 r wx 2 /a) . (3) 

Let (t,r) — > (t + r, f + R) be a uniform translation in space and time, where r is an 
arbitrary time and R is a lattice vector. The classical equations of motion remain invariant 
under these transformations; whereas the Schrodinger equation does not, the reason being 
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the space and time dependence of the gauge potentials. Nevertheless, quantum dynamics 
of the system remain invariant under the combined action of space-time translations and 
gauge transformations. Following Ashby and Miller J7| we define the electric and magnetic 
translation operators 

T (r) = exp (-irO ) , 7} (a) = exp (ia Oj) , (4) 

where j = 1,2 and the symmetry generators are written as new covariant derivatives O^ = 
V\i + ''Vj with the components of the gauge potentials A^ given by 

A = (3f-E, 

A = (" + 1/2) Bx 2 - (0 - 1) E x t , 

A 2 = (a-l/2)Bx 1 -((3-l)E 2 t. (5) 

The symmetry operators in Eq. (Q) commute with the operator S in Eq. (|IJ). The electric- 
magnetic operators given by Ashby and Miller include simultaneous space and time trans- 
lations; we deemed it more convenient to separate the effect of the time evolution generated 
by the T to that of the space translations generated by Tj. The following commutators can 
be worked out with the previous expressions 

ho, TTj] = -iEj , [tti, tt 2 ] = -%B , 

[0 , 0,]=^-, [0 1 ,0 2 ]=iB, (6) 

with all other commutators being zero. The commutators in the second line of Eq. (^|) 
are part of the more general Lie algebra of the magnetic-electric Euclidean two dimensional 
group |J. Schrodinger's equation and the symmetry operators are expressed in terms of 
covariant derivatives 7r M and O^, respectively. A dual situation in which the roles of 7r M and 
O^ are interchanged, could be considered. According to Eqs. (^JjlH), the dual problem 
corresponds to a simultaneous reverse in the directions of B and E. 

The symmetry operators in Eq. commute with S but they do not commute with 
each other. We follow a 3 step method to find a set of simultaneously commuting symmetry 
operators. (1) First we consider a frame rotated at angle 8, with axis along the longitudinal 
and transverse direction relative to the electric field. An orthonormal basis for this frame is 
given by = (cos9,sin9) and = (—sin9,cos9). We assume a particular orientation of 
the electric field, for which the following condition holds 

p = tan£ = - = -, (7) 

where m\ and m 2 are relatively prime integers. This condition insures that spatial periodicity 
is also found both along the transverse and the longitudinal directions. Hence, we define 



a rotated lattice spanned by vectors bi = bei and br = bir where b = ay ml + m\. The 

spatial components of the symmetry generator O are projected along the longitudinal and 
transverse directions : Ol = t^-O and Ot — er-O. It is readily verified that [O , Ot] = 0. 

(2) For the rotated lattice, we regard the number of flux quanta per unit cell to be a 
rational number p/q, that is 
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# 3 **=E. (8) 

We can then define the extended supperlattice. A rectangle made of q adjacent lattice cells 
of side b contains an integer number of flux quanta. The basis vectors of the superlattice 
are chosen as qb^ and bj<- Under these conditions the longitudinal and transverse magnetic 
translations Ti(qb) = exp(iqbOi) and T<r(b) = exp (ibOx) define commuting symmetries 
under displacements qbi and bx- 

(3) We observe that T and T L (qb) commute with T T . Yet they fail to commute with 
each other: T (r) T^qb) = Ti(g6)T (r) exp (—iqbrE). However the operators T and T^qb) 
will commute with one another by restricting time, in the evolution operator, to discrete 
values with period 

2vr 1 / 6 \ 
qbE p\v d J 

where = E/B and we utilized Eq. flD to write the second equality, b/vj, is the period of 
time it takes an electron with drift velocity to travel between lattice points. Hence, the 
meaning of (^j is that the ratio of the Stark ladder spacing (bE) to the Brillouin zone for 
the quasienergy (27TVd/b) is given by the rational number <fi = p/q. 

We henceforth consider that the three conditions (|7|), (H), and (|9|) hold simultaneously. 
In this case the three EMB operators: the electric evolution % = Tq(tq) and the magnetic 
translations = Ti,(qb), and 7^ = T^ib) form a set of mutually commuting symmetry 
operators. In addition to the symmetry operators, we can define the energy translation 
operator M 

T E = exp (- A) > (10) 

that produces a finite translation in energy by 27t/to = qbE. Te commutes with the three 
symmetry operators but not with S. Its eigenfunction exp (—iqbEC,) defines a quasitime £ 
modulo r . 

Having defined %, and 7^ that commute with each other and that also commute with 
S, we can look for solutions of the Schrodinger equation characterized by the quasienergy £ 
and quasimomentum k^ and kx quantum numbers according to 

To \& = exp (—ir£) ip , 
T L ^ = exp (ikj^qb) ip , 

T T * = exp (ihj-b) ip . (11) 

In particular, the previous relations imply that if a simple space-time translation acts on 
the wave function, the latter satisfy the generalized Bloch conditions 

*(t + r , f) = exp{-zr {£ + A (t, f))} *(t, f) , 
(t, r + qb L ) = exp{iqb (k L + A L (t, f))} *(t, r) , 

(t, f + frr) = exp{i6 (fcr + A T (t, r))} (t, f) . (12) 
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We now find convenient to apply a transformation to new variables given by 



Xq — — t , Pq — Oo , 
E 

Xi = ir T + m— , Pi = 7T L / B , 

X 2 = O l /B, p 2 = o T + m^, (13) 

that satisfy [X^, P M ] = i for /i = 0, 1, 2. The explicit relation between (X^ , P M ) and (x M , p^) 
can be worked out using the definition of n^, and and Eqs. (0) and (|5|). The transforma- 
tion is not canonical because of commutator [Po, X 2 ] = iE/B, all others being zero. Applied 
to Eq.(|l]) the transformation yields for the Schrodinger equation 



B 2 Pf + X\ mv d 



M>, (14) 



2m 

where x\ and x 2 have to be written in terms of the new variables as 

For [/ = the dynamics is cyclic in coordinates X 2 and Xq. The appearance of the time 
variable Xo introduced by the periodic potential suggests that an approximation is required 
to solve the problem (e.g. adiabatic approximation). However, as we show below the use 
of an appropriate representation makes such an approximation unnecessary. We adopt 
the P , P\, P 2 representation ^e,k L ,k T (P) = (Poi P\-> Pz\^, &t) with quasienergy £ and 
quasimomentum and k^. In this representation the X^ operators act as X = id/dP^, 
Xi = id/dPi and X 2 = —i(E/B)d/dP + id/dP 2 . It is readily clear that the substitution 
of these relations in Eq. (|1~5|) eliminates the d/dP$ contribution. For finding solutions of 
Eq.(|TJp we split the phase space (X^, P M ) in the (X%, Pi) and the (X , P ; X 2 , P 2 ) variables. 
For the first set of variables we choose a set of basis functions in Pi given by the Landau 
wave functions Xn(Pi); these functions yield exact eigenvalues of the kinetic part of the 
right hand side of Eq. (|T4"D : (n + l/2)u c , with the cyclotron frequency u c = B/m. For the 
subspace generated by the variables (X , P ; X 2 , P 2 ), we notice that the four operators (Tq, 
T£, Trfi, T^) with all possible integer values of (i,j,k,l) form a complete set of operators. 



The demonstration follow similar steps as those presented by Zak in reference [|1(| • Hence a 
complete set of functions, for the subspace (X , P ; X 2 , P 2 ), is provided by the eigenfunctions 
of the operators (T l , 7^, T^, 7^), we write them down and verify their correctness: 

<P kL , kT , e ,t(P , Pa) = E c m e^ bE ^-^e l ^ mkL 5 (P - S - IqbE) 5 (p 2 - k T + ^-m) , (16) 

l,m 

where we defined a = 1/0 = q/p. It is easy to check that this function automatically satisfies 
the first and third EMB tranlations in Eqs. (|TTD, whereas the second equation is satisfied 
by imposing the periodicity condition c m+p = c m . In addition <p is also eigenfuction of the 
energy translation operator Eq. (|T0|) with eigenvalue exp(—iqbEC,). The wave function \l/ 
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is then expanded in terms of \ n and <fik L ,kT,e,£- The operator Te is not a symmetry of the 
problem, so we have to multiply by a coefficient and add over all possible values of £; 
the resulting wave function can be recast as 

^k L ,k T AP) = E «nXn(Pi) h m _ pl c m e ^ b ^ 5{P -£- IqbE) 5 (p 2 - k T + ^m) , (17) 

n,l,m 

where h m _ p i = d^e l ^ abE ^ pl ~ m \ Of particular interest are the following Bloch conditions 
obeyed by ^ with respect to the eigenvalues: 

* {£, k L + bB, k T ) = * {£, k L , k T ) , 
(£, k L , k T + qbB) = e iqbkL (S + qbE, k L , k T ) . (18) 

These conditions are quite different from those satisfied by the usual Bloch and magnetic 
Bloch functions |J; the second one is not periodic, because in addition to the Bloch phase 
e iqbk L — > forp -\- qb shift leads to the change in energy £—*■£ + qbE. 

Based on wave function fllTD it is possible to work out a complete solution of the EMB 



problem, similar to what has been achieved for the magnetic-Bloch problem [11]. This will 



be analyzed elsewhere |12|]. Nevertheless, the approximated solution wherein the coupling 
between different Landau levels is neglected is interesting enough, and probably more illu- 
minating. Within this approximation the Landau number n is also a conserved quantum 
number and the substitution of Eq. (|17D in (0) yields the secular equation 

Eb. 



^ „ ) Cm fn(& 

2n I 



TL \ ( p-im^m- | imiS m ~ \ 

I TT \* /V"12£ m ~ , -im2S m ~ 

+ U2 a le c m+mi -+- e c m _ n 



(19) 



where we defined: c m = e^ 27r// ^ mfcL /i m c m , S m = (27rm — /c T 6) a, A = exp{— i7rm 1 m 2 cr}, 
/n(c) = e~ na / 2 L n (tto") and L n are the Laguerre polynomials. The quasienergy £ is related 
to the eigenvalue A according to £ — (n + 1/2) oj c + + A. 

The dynamics of the system is then described by Eq. (|19|), a finite difference equation 
with distant neighbor couplings mi and rri2, which also includes a linear term proportional 
to E. This equation generalizes Harper's equation to include the effect of an electric field 
of arbitrary intensity. If the electric field is switched off, it can be set ni\ = 1, m 2 = in 
which case Eq. (|T9| ) reduces to Harper's equation. 

We present results for particular cases when the electric field is aligned along the axis of 
the original lattice (i.e. vn\ = 1, m 2 = 0). It is known that the experimental observation of 
the butterfly spectrum could possible be achieved in lateral surface super lattices QJHJ , given 
that a value for a ~ 1 can be obtained for feasible magnetic fields, due to the large dimensions 
of the unit cell. Hence we select the following values: a ~ 100 nm, Uq = U\ = U2 = 0.5meV, 



m = 0.07m e and E = 0.05 V / cm, that can be satisfied in current experiments [13|. For these 
values, a = 1 corresponds to a magnetic field of 10 T, hence U /uj c ~ 0.05 and the condition 
required for weak periodic potential is satisfied. Figs. (1) and (2) show plots for the scaled 
A spectrum as function of a. We recall that for E = the spectrum for A/ [Uof n (a)} , is 
invariant under the substitution a — > a + N with N an integer. 
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Fig. 1 Quasienergy spectrum for the lowest Landau level. The energy axis is rescaled 
to A/ (Uofo(a)), the parameters selected are: a ~ 100 nm, U = 0.5meV, m = 0.07m e7 
E = 0.05 V/cm and k T = 0. 
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Fig. 2 Quasienergy spectrum of the second Landau level (n = 2). In this case the energy 
axis is rescaled to A/ (Uq), the other parameters are the same as in Fig. 1. 
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In Fig. (1) for the lowest Landau level we observe that for a strong magnetic field the 
Hofstadter butterfly is clearly shown for a in the interval [0, 1]. However, as the magnetic 
field decreases (a increases), some of the fine-grained structures of the spectrum are smeared 
by the effect of the electric field. A weak distorted replica of the butterfly is still observed 
for cr in the interval [1,2]. Finally, for bigger values of a, the spectrum is strongly distorted 
by the electric field and the butterfly "flies away". At this limit the effect of the periodic 
potential is negligible and the spectrum consists of equally spaced levels with separation 
crEb, i.e. a Stark ladder. For the n = 2 Landau level (Fig. 2) three smeared replicas of 
the butterfly can be observed, the narrow waist segments of the plot corresponds to the flat 
band condition given by a = 7 n /7r, with 7 n a zero of the Laguerre polynomial. 

In conclusion, the electric and magnetic translation symmetries are utilized to analyse 
the EMB problem. Bloch functions are derived and their properties established in Eqs. fll2"|) 
and (0). The system is governed by Eq. (|19|), the spectrum of which interpolates between a 
butterfly-like structure and a Stark ladder structure. This equation offers a very interesting 
model, susceptible of analysis in terms of dynamical systems. We finally remark that the 
present formalism should set the basis for the study of Hall conductivity beyond the linear 
response approximation. 
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